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SUMMARY 


A FORTRAN subroutine Is described which provides the capability to plot 
three-dimensional data on an interactive cathode-ray-tube computer terminal i 
conventional plotter. The plotted data, which must be lescrlbed In terms o' 
a series ot two-dimensional turves, appears to form a surface In three- 
dimensional space. Features ot the subroutine include a capability for hutdei 
line computations. User instructions and sample programs are described. 

INTRODUCTION 

The subroutine Ot PTHVW was developed In connection with a need for plotti 
wind-shear profiles. In which the magnitude and direction o f the wir\! varied 
with both altitude and distance from the runway. Convention^ two 11mens1on.il 
plots, which only display function's of one Independent var able, were not 
adequate for presenting the spatial distribution of the wind velocities. The 
P1PTHVW three-dimensional plots are Intended primarily for obtaining a qua 1 
Itatlve anpralsal of the data rather than quantitative vd if" . Axes and sea It 
are thus omitted from the p'uts In favor of picture cla 1- t 

The general- pur pose design of the subroutine makes It suitable for plotting 
any data which can be expressed as a function of two Independent variable- . 
DfPTHVW Is written In FORTRAN and ir 's the set of plot - < b -ou' Mies contained 
the TEKTRONIX Terminal Control System software package (ref. . The user must 
therefore have access to the TIKTRONIX package, although no familiarity with its 


contents Is neress ;ry. The user need only specify in Ms program a set of 
parameters which determines various plotting format options and identifier the 
data to be plotted. 

SYMBOLS AND ABBREVIATIONS 
cathode-ray tube 

name of the three-dimensional plotting subroutine 
a scientific programming language 
Langley Research Center 
Network Operating System 

PROGRAM THEORY AND DESCRIPTION 

A function of twc variables can be represented by a surface in 
three-dimensional space. The two independent variables define a point on a 
base plane, and the value of the function is the height of the surface above 
'or below the plane at that point. If one of the Independent variables is fixed, 
and the other is permitted to assume a range of values, the function will 
represent a "slice" of the three-dimensional surface. We are then restricted 
to a curve on the surface, and the curve is in a plane, i.e., two-dimensional. 

If the process is repeated by incrementing the value of the fixed variable, the 
entire surface may be reduced to a set of curves that can be plotted. 

The subroutine DEPTHVW depicts three-dimensional in1r>miation which is 
described in terms of these surface slices. The slices are plotted slightly 
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offset from each other to give the appearance of a surface (fig. 1). Using 
a greater number of slices, which brings the surface lines closer together, 
generally Improves the surface-like appearance. A representation of the base 
plane and a zero-reference line are plotted to aid the viewer in visualizing 
the date surface (fly. 2). 

The user must provide his own data slices (In the form of a FORTRAN 
data array) one at a time to the subroutine DFPTHVW, along with several plot- 
control parameters. The subroutine then plots the data at the user's cathode- 
ray-tube (CRT) terminal, or alternatively generates a vector file which may be 
used for off-line plotting. User specification of the control parameters 

determines the disposition of hidden lines, endpoint connections , and various 
other options, each of which Is explained in the next section. 

CRT terminals which are compatible for use with the subroutine Include 
the TEKTRONIX 4010, 4012/13, and 4014/15 Computer Display Terminals. NASA 
Langley Research Center (LRC) users may also convert the subroutine output to 
a format suitable for off-line plotting by following the procedures outlined 
in appendi x A. 

PROGRAM USE 

Prospective users must have access to the DEPTHVW subroutine and the 
TEKTRONIX Terminal Control System software package. For Langley Research 
Center users, the following job control statements are required (ref. 2 and J): 
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GET (DVWt IB/UN • 474795N) 

ATTACH (LIBFTEK/UN - LIBRARY, NA) 

LDSET (LIB • OVWLIB) 

LDSET (LIB ■ LIBFTEK) 

DEPTHVW and two supplementary subroutines (a hidden-line routine and an 
endpoint -connection routine) are contained In th'- library file DVWl IB, while 
the TEKTRONIX routines are contained in the library file l IBFTfK. The core 
storage requirement for DFPTHVW and the necessary T1KTRONIX modules is 
approximately 65k octal words. 

For users who do not have access to OVWLIB, a FORTRAN source listing of 
the DEPTHVW routines is provided In appendix B. Langley Research Center users 
may also access tne DEPTHVW source statements by using the following conwand: 

GET (DVSORCE/UN • 474795N) 

The user's FORTRAN main program must assign values to several plot 
‘control parameters, create or access the data arrays, and perform a 
subroutine call to DEPTHVW. The syntax of the subroutine call Is: 

CALL DEPTHVW ( VL F .MODE ,X,Y,M.N,XMIN, XMAX , YM I N , YMAX , I H I DF , I BAUD ) 

Each of the subroutine arguments has a specified type, purpose, values 
and limitations. For all arguments except X and Y, values must be specified 
prior to the first call to DFPTHVW and should not be altered on subsequent 









The use of each argument Is explained below. Also appendix C contains sample 
programs which were used to plot figures 1-10 on a TEKTRONIX 4014 Computer 
Display Terminal . 

VLF. - Viewpoint Location Factor Is a real variable which may take on 
values between 0.0 and 1.0. The purpose of this parameter Is to vary the 
amount cf overlap of slices In the vertical direction, t*us altering the 
appa-ent viewpoint. Low values for VI F cause large overlap and the viewing 
direction appears to be almost parallel to the base plane (flq. 2). Values 
near 1.0 cause liitle overlap and the viewer appears to be looking down on top 
of the data surface (fig. 3). The value of VLF required to get the best 
viewpoint will v«ry with the number of curves, the smoothness of the data, and 
user preference. In general, the greater the number of curves, the smaller 
VLF must be to achieve the same viewpoint. 

MODE . - This Is an Integer variable with possible values 0, 1, 2 and 3. 
This parameter specifies the method of handling endpoints for each data slice. 
A value of 0 specifies no endpoint connection (fig. 4), 1 specifies connection 
to the base plane (fig. 5), 2 causes endpoints to be connected to each other 
and to the base plane (fig. 6), and 3 causes endpoints to be connected only to 
each other (fig. 7). 

X, Y.- These are real arrays containing the data for one slice. Data must 
be sent to the subroutine one slice at a time in order to reduce subroutine 
storage requirements. Thus if twenty curves are to be plotted, the subroutine 
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must ho called twenty times with the X and Y arrays containing data for a new 
slice each time. The x-coordlnates are In the X array and the corresponding 
y-coordlnates In the Y array. Thus X and Y are of equal length. The values 
In the X array should he In strictly Increasing order If the user Intends to 
use the hidden-line capability. Also, the values In the Y array should he 
within the minimum and maximum values specified In YMIN, YMAX or clipping will 
result (especially on the first and last curves), The values In the X array 
must he within the limits specified In XM1N, XMAX. On the terminal screen, 
XMAX Is In the horizontal direction and YMAX Is vertical. 

M. - This Is an integer variable which specifies the length of X and Y 
arrays, l.e., number of data points. 

N. - This Is an Integer variable which specifies the number of curves 
(slices) to he plotted. 

XMIN, XMAX. - These are real variables describing the minimum and maximum 
values of the x-coordlnate scale. Data In the X array should be within these 
limits. Values of XMIN or XMAX mav be negative, zero, or positive, but XMAX 
must be greater than XMIN. The values for XMIN and XMAX should be chosen as 
close to the actual data limits as possible. If the data exceeds these values 
the curve will extend off the base plane, and other unpredictable results may 
occur. On the other hand, If the range of XMIN and XMAX far exceeds the range 
of data values, the screen area Is not utilized efficiently. 

YMIN, YMAX. - These are real variables describing the minimum and maximum 
values of the y-coordlnate scale. Values may be any real number within 
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hardware limitations, but YMAX must be greater than YMIN. Clipping will 
result, especially In the first and last curves due to screen height limita- 
tions, If the data In the Y array exceeds these limits, 

IHIDE. - This Is an Integer variable which specifies the method of 
handling hidden lines. Possible values of this variable are 0, 1 or 2. A 
value of 0 specifies no hidden-line computation Is desired, l.e., lines which 
would normally be hidden will be plotted as solid lines (fig. 8). A value of 
1 causes hidden lines to be plotted as dashed lines (fig. 9). A value of 2 
causes hidden lines to be Invisible (fig, 10). 

If the hidden-line capability Is used (that is, IHIDE nonzero) the user 
should have at least 100 data points per curve. Fewer points than this will 
reduce the effectiveness of the hidden-line checking algorithm. For example, 
a supposedly invisible line will remain visible for some distance even after 
going behind a surface. Also, the X array should be arranged in strictly 
Increasing order. Finally, the hidden-line routine checks only the previous 
10 curves; if a line Is hidden only by a slice 11 or more curves back, it will 
be visible. Despite these limitations, the hidden-line capability is useful 
In reducing the confusion of several curves close together on the screen. 

IBAUD_^ - This Integer variable specifies the computer terminal trans- 
mission rate In characters per second. 
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CONCLUDING REMARKS 


DEPTHVW Is a subroutine which plots a series of data points so that they 
appear as a surface In three-dimensional space. It does not draw a continuous 
surface, but Instead plots a series of two-dimensional curves spaced properly 
so as to give the appearance of a surface. 

The subroutine permits a qualitative, easily Interpreted visualization 
of data that depends on two Independent variables. It has several limitations, 
however. It was not designed to depict complex three-dimensional objects such 
as structures. Also, to nuke full use of Its features, the data must be In 
strict order and within specified limits. 

The subroutine !s still In a developmental stage, and may contain 
deficiencies that will show up with the use of different data. The authors 
solicit user consents on the subroutine In order that further Improvements 
can be made. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 1 , 1978 
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APPENDIX A 


USING DEPTHVW WITH LRC GRAPHICS POSTPROCESSOR 

DEPTHVW may he used to create three-dimensional plots without an 
Interactive CRT terminal, utilizing the IRC graphics postprocessor to creat** 
Varian or fa 1 Comp plots, The following steps are necessary! 

A. Changes to Main Program 

Access to the LIBFTEK co.;mon block /JTB/ Is required. The user must 
Include the following statements In his main program: 

COMMON/JTB/ NFR, JREQ, DUMMY (8) 

JREQ » 1 

The first Is a non-executable statement and must occur before any execuLalle 
statements In the program. 

Setting the value JREQ * 1 In /JTB/ causes the plot vector file to be 
written to file TTPE16 Instead of file OUTPUT (the CRT screen) In order that 
the data may be plotted off-line. In this case, the plot will not appear on 
the screen. If the program Is being run at a CRT. The user may also obtain a 
plot on the screen by using the following batch control statements: 

ASCII. 

LNH.F - T TPE 16. 

R 


B. Plot Vector File Conversion 


The plot vector file must be converted from TEKTRONIX CRT format to a 
format acceptable to the LRC graphics postprocessor. This Is done with 
procedure file MCOMP, and can be done Interactively or In batch. Control 
statements necessary are (TTPE16 Is a local file)* 

GET, MCOMP/UN * 075768N 

REWIND, TTPE16, MCOMP. 

CALL (MCOMP) 

REWIND, SAVPLT, 

MCOMP converts TTPE16 to a LRC graphics vector file which Is stored In SAVPLT. 

C. Postprocessing 

NOS does not allow postprocessing to be done Interactively. The vector 
file conversion described above can be done Interactively and postprocessing 
In a separate batch job, or the two steps can be combined In a single batch 
job. 

Note that If conversion and postprocessing are done *s a single batch job, 
the last two control statements In procedure file r „jMP should be deleted 
before It Is called. That Is, the statements, 

EXIT. 

REDUCE. 
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should be deleted from the user's copy of MCOMP before It Is used In the 
combined Job described above. 

Postprocessing requires the following batch control statements (SAVPLT Is 
a local f i le) i 

ATTACH (PLOT/UN ■ LIBRARY) 

PLOT. "Device Name". . , 

Reference 4 gives details on the IRC plotting devices and how to control them 
with the PLOT control card. 


11 


APIUNDJX n 

SOURCE LISTING OF DEPTMVW SUBROUTINES 


SUBROUTINE DEPTH VW < VL^ 1 • NMOPE 1 t XP , VP , M 1 , N1 , 
•XMIN1 tXMAXl # VMlNl tVMAXl t IHIDE1 t IBAUD) 
DIMENSION XP(Ml) t YP(Mn 
COmmON/PARAm/ vlf«nmooe« imioe* 
•IRLNK,TESTfMFXC.nELX.OELV f M,N. It J1 « 

•IBLNK? f IFLGtXMAX 
DATA IflLNK/O/t IE IPST/1/ 

C 

lEdFIBST.NE.n on TO 100 
C LOAD THE COMMON BLOCK 

VLF * VI.F 1 

NMODE*NMODEl 

IHIOF^IHIOEI 

M*M 1 
N«N1 

XMA X ■ XHA V 1 
XMINaXMINl 
VMIN«YMIN1 

ymax*ymax 1 

C INITIALIZE ANO CALCULATE winoo* size 

C', INI TT ( IMAun) 

. M 

iWOTHslO ? A • ? / 3 

VHGHTb (VMAX-YMIN) »N»VLF* (YMXX-YMIN) 

VWDTH* 1 ,B»(XMAX-XMIN) 

vbot*ymin 

VLET*xmin 

CALL VKIN00( VLFT tVWOTHt VBOT, VHGHT) 

C OR Aw BASE PLANE 

IHGHT« (VH«HT-( VMAX-YMIN) )/VHGHT»7B0 
CALL MOVABS(O.O) 

CALL ORWARSC IWOTH.O) 

CALL DPWABS ( 10?*t IHGHT) 

IFIPST=0 

100 IF(l.GT.N) GO TO 999 

HFAC«(XP(M)-XP(1 1 ) /100. 

TEST*XP ( 1 ) 

CALL MOVEA (VLFTt VBOT) 

C DETERMINE offset for this curve 

0ELXSVW0TH/3.* ( 1-1 ) /(N-l ) 

DFLT= ( VHGHT- ( YMAX-YMIN) 5 • < I - 1 > / < N- 1 ) 

Ir LG*0 

00 ZOO J=1.m 

C DETERMINE COORDINATES ANO GRAPH POINT 

XCPDsXPI J) *DEL X 
YCRD«YP( J) « DEL Y 
IF(J.EO.l) XCD?»XMIN*DELX 
IF(J.EO.l) YCDpsYCRO 


UKK.UNAL 

OF POOR 


t Avi Li R’ 

u UAAU^ 


n o on 


ORIGINAT, PACK I. 
OF POOP QUALITY 

IF(J.LQ.l) CALL MOVEA (XCRDtYCRD) 

Jl*J 

IF(IHIDE.NE.O.A.XPU) .GE.TEST) CALL HlDOEN(XCRD*YCRD) 
CALL D-SHA<XCPD.YCPD*I£5LNK) 

200 CONTINUE 
C ENOPOINT HANDLING 

DOWNl*YP( JU-VMIN 
CALL MOVEA (XMAr.^OELX.YCPO) 

IF(I.EQ.l) OOWN?«YP(l)-YMlN 

CALL ENDPT (XCD1 tYCOl » DOWN 1 * DOWN2 *XCD?*YCO?) 
XC01*XMAX*DFLX 
YC01«YCRD 
I * I ♦ I 

999 RETURN 
END 

SU3R0UT I NE ENDPT < XCD 1 # YCD1 .DOWNl ♦DOWN?, 

•XCD?* YCD2) 

COMMON/P AM AM/ VLF »NMODE« IHIDE* I BLNK « TEST « HF AC • OFL X » 
•OELY.M.N. I , J, IHLNK2, IFLGtXMAX 
DATA IBLNK3/0/ 

MODE 0 IS DO NOTHING 
IFLG«1 

IF (NMODE .FO.O) GO TO 99 

MODE 3 IS CONNECT ONLY FIRST ANO LAST CURVES 
TO 8 ASF I CONNECT ENDS TO EACH OTHER 
IF(NMODE.EO,3.A.I.NE.1.A.I.nE,N) GO TO ?0 
C MODE 1 IS CONNECT FNDS TO RASE ONLY 

CALL DRAwR(O.t-DOwNl) 

CALL MOVER ( 0 • * DOWN 1 ) 

IF(NMOr*E.FO.n 60 TO 99 
20 IF(I.EO.l) GO TO 30 
C CONNECT FNDS TO EACH OTHER 

CALL DPAWA (XCD1 t YCD1) 

CALL MOVEA (DELX*XMAX tOELY) 

CALL DRAWA (XC04, YCD4) 

30 CALL MOVEA (XC02.YCD2) 

CALL HIDDEN(XCD?.YC02) 

I8LNK?=IBLNK 
IF(I.EQ.l) GO TO 50 

C CHECK FOR HIDDEN END-LINES ON LEFT SIDE 

40 IF(IBLNK2.NE O.O.IBLNK3.NE.O) GO TO 70 
GO TO 80 

C DRAW CONNECT TO BASE ON LEFT SIDE IF FIRST CURVE 

50 CALL DPAWH (0. t-DOWN2) 

CALL MOVER (0. .D0WN2) 

C PUT PRESENT ENDPT COORD. IN OLO COORD. LOCATIONS 

60 XCD3=XCD2 
YCD3= YCD2 
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1HLNK 1»1«l NK? 
acp4*de l a ♦ ama a 
vcP4*of iv 
00 TO RR 

C OFTFRmINF MOW MUCH OF LEFT FNO-LtNF IS MlDOFN 

C AND DRAW FNO-LtNF CORRECTLY 

7 0 IF(IBlNK?.NF.n.A.IHLNK3.NF,0> 00 TO AO 
XRT-XCO? 

XlT*XCD3 

V TOP* VCD? 

Y ROT ■ Vf D 1 
I F L 0 ■ 1 

DO 7R ITST*1.10 

AMID* (KPT^ALT)/?. 

YMin*(VTOP*vHon/?, 

CALL MIDI>FN( AMID.YM1D) 

IF lim NK ,f (1,01 00 TO 75 
IF < IHLN5?.FO.O) 00 TO 76 

74 XRT*XMIP 
VT0P*VMIP 
00 TO 7R 

75 IF (Ifll NK?.FQ.O) 00 TO 7* 

76 ALT*AMID 
YBOT *VMI D 

7R CONTINUE 

CALL DASHA ( amid. YM1D* IHLNK?) 

HO CALL DASHA (ACD.T#YC03* IBLNKH 
00 TO 60 
RR RETURN 
FND 

SUBROUTINE HIDDF N( ACRD# YCMD) 

D IMF NS I ON ACD ( 10,tl0»#VCP(10.110) 

COMMON /PAR AM/ VL F » NMOOF , I M I DE • I HL NK • TF 5 T t MF AC t DFl A, 
•OFLY.M.N. I ♦ J, I RINK?. I FLO. AMAA 

INITIALIZE IF FIRST CURVE. FIRST POINT 
IFUI.NF.l.OR.J.Nf.n.A.IHIOE.NE.O) 00 TO S 
00 4 JHL A * 1 # I 0 
DO 4 KHL A«t ,110 
ACPI JBLA.KBLA) *-1000000, 

YCP( JHL A, KHL A I *-1000000, 

4 CONTINUF 

5 IFIJ.EO.I) OOT 030 

C DFTFRMINF IF POINT IS HIOOFN 

6 IF ( I . F 0 , I ) 00 TO ?0 
DO 16 L *? • 1 0 

IF(XCR0.GT.XC0IL»IST)) 00 TO 16 
IF < ACRD.LT ,ACO<L. ’ M 00 TO 16 
00 10 K * 1 # L S T 
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OK I(J INAL PACK IS 

OK POOR OIJAU1T 


IFIXCRD.L! .XCO<L«K)> GO TO IS 
10 CONTINUE 

15 IF (YCR0.LT.YCD(L*K) ) 60 TO ?\ 

16 CONTINUE 

C SET MIODCN FLAG 

?0 IRLNK«0 
GO TO 22 

21 IBLNK*-| 

ifuhidf.fq.p I RLNK * 3* 

22 IFHFLG.NF.fl) GO 10 99 
XCDU* ICNTjnXCRI) 

YCO ( 1 . ICNT) *VCRO 
IF(J.NE.l) rFSTsTFST^HTAC 
ICNT * I CNT ♦ 1 
GO TO 99 

•JO IF(IFLG.NE.O) GO TO 6 
IF(I.NE.l) LST*ICNT-1 
ICNT * 1 

IF ( IFLG.NF.O) GO TO 6 

C LO AO PPESFNT COORD. INTO OLD COORD. LOCATIONS 

DO 35 KOMtllO 
00 35 L* 1 • 9 
LO* 1 1 -L 

XCO(LO.KO) =XCD (LO- 1 # KO) 

35 YCD(LO.KO) *VC0(L0-1 »K0) 

GO TO 6 
99 RFTURN 
END 
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APPENDIX C 

SOURCE LISTING OF SAMPLE PROGRAMS 


NOTE: Program CHECK1 was used to qenerate figures 1, 8, 9 and 10. 

Program HYPLOID was used to generate figures 1 - 1 . 

. i 

PROGRAM CHECK1 < INPUT. OUTPUT) 

C THIS PROGRAM PLOTS A SPIKE 

DIMENSION X(1001)«Y(1001) 

PRINT* ♦"NMOOF«VLf ♦ I HIDE" 

REAO*.NMOOE.VLF. IHIDE 
C CREATE DATA ARRAYS 

DO 20 1*1.21 

Z*EXP (FLOAT C (1-1 1 •(21-1) > - 1 00 > ) 

DO 10 J»1.1001 
X ( J) *FLOAT ( J-l ) /100. 

Y( J)*Z*EXP (X ( J) • ( 10.00-X ( Jl n«75./EXPf5.**?) 

10 CONTINUE l 

TRANSFER TO DEPTHVW WITH DATA ARRAYS AND CONTROL 

PARAMETERS 21 TIMES 

CALL DEPTH VW (VLF.NMOOE.X.Y. 1001* 21 t 
•0. . 10. t- 10. .SO, . I HIDE* 120) 

20 CONTINUE 

WAIT UNTIL USER PRESSES ‘'RETURN' 1 
CALL TPAUSE 
END 


PROGRAM HYPLOID < INPUT. OUTPUT) 

THIS PROGRAM PLOTS A HYPERBOLIC SURFACF 
DIMENSION X ( 1000) ,Y( 1000) 

DATA XMIN.XMAX.YMIN.YMAX.M.N/-100..100..-10000,. 

• 10000 .* 1000 . 20 / 

PRINT". "MODE.VLF. IHIDE ?" 

RE AO*. MODE .VLF. I HIDE 

DATA CREATION AND CALL TO DFPTHVW - - N CURVES. 

M DATA POINTS PER CURVE. N CALLS TO PLOTTFR 
DO 20 1*1. N 

FIX=100.-(I-1)*200./1R. 

DO 10 J = 1 • M 

XCJ)=-10C-M J-l ) *200./R9R. 

Y(J) = (X( J)**2-FT '.**2) 

10 CONTINUE 

CALL DERTHVW ( VLF .MODE. X. Y.M.N, XMIN. XMAX. 
•YMIN.YMAX. IHIDE. 120) 

20 CONTINUE 
C WAIT UNTIL USER PRESSES ''RETURN'* KEY 

CALL TPAUSE 
END 
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